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We perform state-of-the-art large-scale shell-model calculations of the structure factors for elastic 
spin-dependent WIMP scattering off ^^^'^^^Xe, ^^^I, ''^Ge, ^^F, ^^Na ,2^A1, and ^^Si. This compre- 
hensive survey covers the non-zero-spin nuclei relevant to direct dark matter detection. We include a 
pedagogical presentation of the formalism necessary to describe elastic and inelastic WIMP-nucleus 
scattering. The valence spaces and nuclear interactions employed have been previously used in 
nuclear structure calculations for these mass regions and yield a good spectroscopic description of 
these isotopes. We use spin-dependent WIMP-nucleus currents based on chiral effective field theory 
(EFT) at the one-body level and including the leading long-range two-body currents due to pion 
exchange, which are predicted in chiral EFT. Results for all structure factors are provided with 
theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei. 
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I. INTRODUCTION 

About 25% of the energy density of our Universe con- 
sists of dark matter, a form of matter that rarely inter- 
acts with baryons and has eluded direct observation so 
far [l|, |2| ■ This large-scale problem is closely connected 
to new physics at the smallest scales, because dark mat- 
ter candidates arise naturally in extensions of the Stan- 
dard Model of particle physics [3] . Prominent dark mat- 
ter candidates are weakly interacting massive particles 
(WIMPs) that are predicted as the lightest supersymmet- 
ric particles (typically neutralinos) in supersymmetric ex- 
tensions of the Standard Model. WIMPs are especially 
promising candidates, because they account naturally for 
the dark matter energy density established by observa- 
tions [2|- Moreover, WIMPs interact with quarks, and 
thus with baryonic matter, opening the door to direct 
dark matter detection via elastic scattering off nuclei [J] . 
Inferring properties of dark matter from direct detection 
therefore requires detailed knowledge of the structure fac- 
tors for WIMP scattering off strongly interacting nuclei. 

In this work, we focus on spin-dependent (SD) WIMP 
scattering 5] , which is relevant because neutralinos carry 
spin. The detection of elastic SD WIMP scattering has 
been the goal of several past and ongoing experiments [6|- 
[l3 |. using different non-zero-spin nuclei as target, but 
so far without evidence. Evaluating the response of nu- 
clei to WIMPs is challenging. First, it requires matching 
the WIMP-quark couplings in a particular supersymmet- 
ric model to WIMP-nucleon currents. Because quantum 
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chromodynamics (QCD) is nonperturbative at low en- 
ergies, this is best achieved using effective theories |13l - 
|l6[, in which spin- independent and SD interactions gen- 
erally enter in leading order. Second, the WIMP-nucleus 
response requires reliable nuclear-structure calculations. 
This is especially important for SD interactions, because 
the response depends on how the spin of the nucleus is 
distributed among nucleons (due to attractive interac- 
tions most of the nucleons pair to spin zero). For the 
isotopes of interest [6l-[l3|. this involves medium- mass 
to heavier nuclei and is a challenging many-body prob- 
lem. Previous calculations of SD WIMP scattering off 
nuclei 0, [l7H23j have relied on phenomenological WIMP- 
nucleon currents, and are based on nuclear-structure cal- 
culations that can be improved with recent advances in 
nuclear interactions and computing capabilities. This 
work presents progress on these fronts. 

The typical momentum transfers involved in WIMP 
scattering off nuclei are low and of order of the pion mass. 
In addition, the typical momenta involved in low-energy 
nuclear structure are similar. At these momentum scales, 
chiral EFT provides a systematic expansion in powers of 
momenta Q for nuclear forces and for the coupling to ex- 
ternal probes, based on the symmetries of QCD |24|.|25|. 
In addition to the coupling through one-body (lb) cur- 
rents, generally at leading order, two-body (2b) currents 
enter at higher order and are quantitatively important. 

In previous work [16] we have derived the currents for 
SD WIMP scattering off nuclei based on chiral EFT, in- 
cluding lb currents and the leading long-range 2b cur- 
rents due to pion exchange, which are predicted in chiral 
EFT. As an application, we focused on the scattering off 
129,131-j^g^ as they provide the most stringent limits for 
WIMP coupling to neutrons [2g. Our results have re- 
cently been adopted as benchmark for the XENONIOO 
SD WIMP-nucleon cross-section limits [I4I . 

More generally, two-body contributions to weak neu- 



tral currents have been shown to be key for providing 
accurate predictions of neutrino-deuteron scattering at 
solar neutrino energies for SNO [23, [23|- Weak neu- 
tral currents based on chiral EFT have been explored for 
light nuclei and neutrino breakup in core-collapse super- 
novae |29l - l32l |. and 2b weak charged currents have been 
shown to provide important contributions to Gamow- 
Teller transitions and double-beta decays of medium- 
mass nuclei |33j . Following our previous work [l6|, 
Refs. [33, [33] have reported simple prescriptions to ap- 
proximately include the effects of chiral 2b currents in 
previous calculations of SD WIMP-nucleus scattering. 

This work expands Ref. ^16,] by presenting state-of-the- 
art large-scale shell-model calculations that describe the 
non-zero-spin states of all isotopes that are experimen- 
tally relevant for SD WIMP direct detection: ^^aasixe, 
12^1, 73cg^ i9p^ 23^^, 27A1, and ^^Sl The nuclear- 
structure calculations are performed with interactions 
and valence spaces that have been tested in these mass 
regions. Based on the calculated ground states, we pre- 
dict the structure factors for elastic SD WIMP scattering, 
including chiral lb and 2b currents at the same level as 
in Ref. [16J. We provide theoretical error bands due to 
the uncertainties of WIMP currents in nuclei. 

The outline of this article is as follows. In Sect. |lll 
we derive the WIMP currents in nuclei based on chi- 
ral EFT. All microscopic inputs needed to compute the 
structure factors of SD WIMP-nucleus scattering are dis- 
cussed in Sect. IIIII Combined with detailed Appendices, 
this includes a pedagogical presentation of the formalism 
necessary to describe elastic and inelastic WIMP-nucleus 
scattering. In Sect. IIVI we present large-scale nuclear- 
structure calculations that describe the nuclei relevant 
for SD WIMP direct detection, and compare our results 
to experiment. We then calculate the structure factors 
for elastic SD WIMP scattering for all cases using chi- 
ral EFT currents. We discuss in detail the role of 2b 
currents and their uncertainties; the contributions of dif- 
ferent multipole operators to the total response; and the 
issue of proton/neutron versus isoscalar/isovector decom- 
positions of the structure factors. Finally, we summarize 
in Sect. |V] and give an outlook for future improvements 
of the nuclear physics of dark matter detection. 



where Gp is the Fermi coupling constant, and J^(r) 
and j^(r) denote the hadronic current and the leptonic 
current of the WIMP, respectively, x is the neutralino 
field, ipq are the fields of q = u,d, s quarks, and Aq 
the neutralino-quark coupling constants. The tempo- 
ral components can be neglected, because the veloci- 
ties of WIMPs are expected to be nonrelativistic with 
v/c ^ 10~^. We also neglect contributions to the La- 
grangian density other than axial- vector currents, such as 
polar- vector currents, which are suppressed by the mo- 
mentum transfer over the nucleon mass p/m [5|. This 
approximation will be studied in a future paper. 

For the WIMP-nucleus response, the SD WIMP inter- 
action couples dominantly to a single nucleon, but also 
to pairs of nucleons. At the one- nucleon level, the quark 
currents are replaced by their expectation value in the 
nucleon, leading to lb axial- vector currents Ji.ib- In the 
nucleus, the currents are summed over all A nucleons: 



] Aqll^qllzi^q > ^ Jj,lb 
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with the isoscalar J^ j^j^ and isovector J^ j^j^ parts. 
The coupling of the isoscalar part is given by [5| 



ao 



{Au + Arf)(Au H- Ad) + 2As/^s , 



(3) 



where Au,Ad,As are defined as V'gTTsV'rj = Aqcr/2, 
with the nucleon spin cr/2. Therefore, J"^^, — ag cr/2, 
and ao receives contributions from the isoscalar combi- 
nation of the u and d quarks to the spin of the nucleon, 
as well as from the s quark. Analogously, the isovector 
coupling can be written as 



ai 



= {Au - Ad){Au - Ad) = {Au - Ad)g/ 



(4) 



where gA is the axial coupling constant. This shows that 
the isovector part J^ -^^^ of the axial- vector WIMP-nucleon 
coupling is identical, up to replacing oi by gA, to the 
axial- vector part of the weak neutral current. 



B. Coupling to one nucleon 



II. WIMP-NUCLEUS INTERACTIONS 

A. Chiral EFT and WIMP currents 

At the WIMP-quark level, the low-momentum-transfer 
Lagrangian density C for SD interactions is taken to be 
an axial- vector-axial- vector coupling [5|, |36[ : 



G 

^^ "72 
G 



£SD ^ ^ / d3rj^(r)J^(r) 



--| / rf^rx775X • Y^ Aq^Pgjj^tPq , (1) 



The weak neutral current was derived within chiral 
EFT for calculations of low-energy electroweak reactions. 
At lowest orders Q° and Q^, there are only lb currents. 
For the isovector part of the axial- vector WIMP-nucleon 
current, this leads to UM 



t3 "'^ „ ^3 



9a{p^ 
9A 



2mgA 



(P • crO P 



(5) 



where t? denotes the isospin, p — Pi — Pi the momen- 
tum transfer from nucleons to neutralinos, and 5a (p^) 
and gp{p^) the axial and pseudo-scalar couplings. The 
momentum transfer dependence of 5a (p^) and gp{p^) is 



due to loop corrections and pion propagators. To order 
Q^, one has (stI 



9A 
9p{p') = 



A2 



-4 



p^ 






(6) 
(7) 



with Aa = 1040 MeV, pion mass m^ = 138.04 MeV, pion 
decay constant F^r = 92.4 MeV, and g^rpn = 13.05. Chiral 
lb currents are similar to the currents used in previous 
calculations of WIMP scattering off nuclei Q . The dif- 
ferences are that the 1/A^ terms were neglected and the 
Goldberger-Treiman relation was implicitly used to write 

^^Jf ' « ^'i\. i - Both present few percent corrections, 
but the former increases with momentum transfer. 

The axial- vector part of the weak neutral currents is 
isovector in the Standard Model, neglecting the strange 
quark contribution to cq in Eq. ([3]). Therefore, higher- 
order Q^ contributions to the isoscalar WIMP-nucleon 
current J^ -^^^ ~ ao <t/2 depend on models of currents in 
the nucleon. To order Q^, these lead to lb currents with a 
form-factor mass-scale ^ Aa |38| and without pion prop- 
agator contributions. Because the isovector 1/A^ terms 
contribute at the few percent level for the typical momen- 
tum transfers in WIMP scattering, we chose to neglect 
higher-order isoscalar current contributions, as opposed 
to introducing a model dependence at this level. 



C. Coupling to t^vo nucleons 



and heavy nuclei, because of phase-space restrictions of 
normal Fermi systems at low energies [40i | . 

The resulting effective lb currents J^^t are derived in 
detail in Appendix[3 We find that the leading long-range 
2b currents lead to three different contributions. First, 
there is a renormalization of the axial coupling jig , 



jf2'b(P.^) 



T? p 



-9Acr^ — J^ 



\{3l2iP,P) 



^^^(^'^))(r^+6l 



Iiip,P)^c, 



(9) 



which depends on the density p and the total momentum 
P = Pi -f p^ (due to the momentum transfer in the ex- 
change term) . Such renormalization was also found con- 
sidering chiral three-nucleon forces as density-dependent 
two-body interactions [41]. Second, there is a contribu- 
tion to the pseudo-scalar coupling. 



J?,^b(p) = -.9A^(p-cr,)p 



2C3 



F^ 4m^ + p2 



(10) 



and third, chiral 2b currents induce a pseudo-scalar-type 
current depending on the total momentum. 



rOff,P 
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{p,P)=-gAfiP■^^)P 

X^/^(p,P)(i(-C4-C3)~-^ 

r^ V3 12?Ti 



(11) 



At order Q^, 2b currents enter in chiral EFT (H]. We 
consider their long-range parts due to pion exchange, 
which are predicted in chiral EFT, and for medium- 
mass nuclei were found to dominate over the short-range 
parts [16i |. Because of their pion-exchange nature, the 
axial- vector part of the weak neutral 2b current is isovec- 
tor, J2b = Ej<j Jfjj with 



•J 12 



9A 



1 



2F^ ml + k^ 

+ 2c3k- (a-lT^ + (T2T^)k 



(c4 + ^)kx (rrx xk)r3 



2m 



k • (o-i - 0-2 )qr^ 



(8) 



where r^ = (n x T2)^, (t^ = {(Ti x 0-2), k = i(p^ -p2 - 
p'l -fpi) andci = i(pi +p'i -P2-P2)- 

As in Ref. [16|, we take into account the normal- 
ordered one-body part of chiral 2b currents. This is 
obtained by summing the second nucleon j over occu- 
pied states in a spin and isospin symmetric reference 
state or core, which we take as a Fermi gas: J^^^ = 
^ (1 — Pij)3fj. The exchange operator P^- includes all 
two-body exchange contributions. Normal-ordering is ex- 
pected to be a very good approximation for medium-mass 



The functions /f(p,P), l2{p,P) and I^{p,P) are given 
by integrals due to the summation over occupied states in 
the exchange term. They can be evaluated analytically, 
and the explicit expressions are given in Appendix [X] 

The contributions from 2b currents in Eqs. (P))- pT|) de- 
pend on the density of the reference state p — 2kp/{3n'^) 
(fcp is the Fermi momentum) and on the low-energy cou- 
plings C3 and C4. For the density p we take the range 
p = 0.10. ..0.12 fm~ , appropriate for the nuclei consid- 
ered (see also Ref. [SJl). The low-energy couplings C3 
and C4 also enter pion-nucleon and nucleon-nucleon in- 
teractions and have been determined from data. Here, 
we consider the 03,04 values from the N'^LO NN poten- 
tials of Ref. dl] (EM) and Ref. ^ (EGM), as weh as 
from the NN partial wave analysis (PWA) [44|. To be 
conservative, we also consider the estimated uncertainty 
in these values expected from higher-order order contri- 
butions, Sc3 = -Sc4 « lGeV~^ [23|. The resulting C3 
and C4 values are given in Tables U and HT] 

In Table HI we study the P dependence of the 2b 
current contribution to the axial coupling, which we 

write as 3l,^^{p,P) = -gA(ri^-^J'^{p,P). We com- 
pare J'^{p,P) for the Fermi gas mean- value P^ = 6fc|/5 
and P = at a density p = 0.10 fm~ and for the dif- 
ferent C3,C4 sets considered. Table |T] shows that the P 
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EGM 
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EGM+5ci 
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PWA 


-4.78 


3.96 


2.77 


2.55 


0.08 


PWA+Sci 


-3.78 


2.96 


2.14 


1.97 


0.08 



TABLE I. Comparison of J" {p, P), which describes the axial 
contribution from the normal-ordered one-body part of the 

long-range 2b currents, J'iJ^iPjP) = -gAcri-^-^J'^(p,P), 

evaluated at the Fermi gas mean- value P = 6fcF/5 and at 
P = for a density p — 0.10 fm~^. The variation is shown 
for all C3,C4 sets considered in this work, and the relative 
variation A J''/ J"^ between the Fermi gas mean- value and P = 
is given. The Ci and J" values are in GeV~^. 
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EM 


-3.2 
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-(0.25.. .0.32) 


0.12.. .0.14 


EM+Sci 
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EGM 


-3.4 
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-(0.19.. .0.23) 
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EGM+5c^ 


-2.4 


2.4 


-(0.13.. .0.17) 


0.09.. .0.10 


PWA 


-4.78 


3.96 


-(0.23.. .0.29) 


0.17.. .0.21 


PWA+5ci 


-3.78 


2.96 


-(0.18.. .0.22) 


0.14.. .0.16 



The 2b current contribution to the pseudo-scalar cou- 
pling, Eq. (jlOp . we similarly write as a momentum- (and 
density-) dependent renormalization Saf{p) [la]: 



Jl5^ip) 



-9A- 



Safip) 



(p • o-'Op . 



with 



Saf{p) = -2c3 



F^ Ami. + P^ 



(15) 



(16) 



The ranges of 5ai and 5a^{p) are given in Table HIl for 
the C3,C4 values and the density range considered. We 
find that 5ai reduce the axial part of the current by 
13 . . . 32%, while Jaf (p) increases the pseudo-scalar part 
of the current by 8 ... 21% at p = m.^-. At lower momen- 
tum transfers this enhancement is weaker, while it is more 
significant for higher p. This is consistent with studies 
of Gamow- Teller transitions and double-beta decays [3j] . 
As discussed in Ref. [1^1 , in addition to the long-range 2b 
pion-exchange currents, there are short-range 2b currents 
for the isoscalar and isovector parts, which are included 
as contact terms in chiral EFT. The isovector short-range 
2b parts only lead to small contributions [3^ . Therefore, 
we neglect short-range 2b currents at this level, which is 
also consistent with neglecting higher-order (short-range) 
lb isoscalar currents, see Sect. IIIBI 



TABLE IL Values for all C3,C4 sets considered of the long- 
range 2b current contributions Sai (axial) and (5aj (p — rTiyr) 
(pseudo-scalar) for the density range p — 0.10...0.12 fm~^. 
The Ci values are in GeV~^. 



dependence is very weak: J'^(p, 0) varies by less than 
14% over the relevant P range. For other densities in the 
range p = 0.10...0.12fm~'^ this variation is even smaller. 
Because 2b currents are a correction to the leading lb 
currents, we therefore set P = in the axial 2b current 
contribution, Eq. ([9]). Similarly, the contribution from 
Eq. (fTTj) is expected be weaker, and we therefore consis- 
tently set P = at this level. 

For P = 0, both /f and I^ lead to ^ 



r{p,P = 0) 



IliP-P-- 
iml 



0) 



^2{P.P 



0) 



= 1- 



Kp 



3m3 / fcp 
3 arctan I 



(12) 



where I{p,P = 0) = 0.58. ..0.60 depends only weakly on 
the density in the range p = 0.10. .. 0.12 fm"^. The 2b 
current contribution to the axial part, Eq. ([SJ, can then 
be written as a density-dependent renormalization 5ai, 



|-Cff,(7 

'i.2h 



-QAtTi-^ 5ai 



with 



5a, = -^rip,P = 0) fi(2c4 - C3) + ^ 
±'4 \o brn 



(13) 



(14) 



D. Combined response 

Combining the lb and the long-range 2b currents to 
order Q^ in chiral EFT (replacing gA by ai for the latter), 
the isovector part of the axial- vector WIMP current at 
the normal-ordered one-body level is given by [16| 



t3 3 

•'iAh+2h ~ 2 '^^'^* 



( 9Aij?) 

V 5A 



5a\ 



9P{P^ 
2mg/ 



Sa^ip) (p-o-,)p 



(17) 



III. 



WIMP-NUCLEUS SCATTERING AND 
STRUCTURE FACTORS 

A. WIMP-nucleus scattering 



The differential cross section for SD WIMP elastic scat- 
tering off a nucleus in the initial state \i) to the final state 
I /) can be obtained from the low- momentum-transfer La- 
grangian density of Eq. ([1]). A detailed derivation is per- 
formed in Appendix |B] The final result is [5| 



da _ 1 

djp ^ (2 J, + l)7rw2 



E EI(/i4°N)l 

Sf,Si Mf,Mi 



^^^ SAip), 



(2 J, + l)t-2 



(18) 



where the sum s/, Si = ±1/2 is over neutrahno spin pro- 
jections, and the sum Mf,Mi is over the projections of 
the total angular momentum of the final and initial states 
Jf,Ji, respectively; v is the WIMP velocity, and Sa{p) 
the axial- vector structure factor. The structure factor 
can be decomposed as a sum over multipoles L with re- 
duced matrix elements of the longitudinal £|^ , transverse 
electric T^'^, and transverse magnetic f^^^^ projections 
of the axial- vector currents; 



L>0 



L>1 



(19) 



The multipole contributions are obtained from the 
WIMP-nucleus currents J'^(r). At the effective one-body 
level, chiral lb and 2b currents lead to (see Appendix |B] 
for the definition of the multipole operators and details 
of the derivation) 

Clip) 



V2L + T 



flo + aiTj 






'^QTrpn^TrP 



<5af(p) 



2mgA{p'^ + ml) 
^L + lMLX + l(pYi) + ^/LML,L-l{pv^)\ , 



(20) 



rrip) = 



V2L + T 



X 



P 



^ - ao -I- airf (^1 ~ 2j^ + 5ai 



x[- 



\^-VLML,L+lipr^) + ^L + lML,L-l{pv{) 



(21) 



rr^'^p) 
^ 1 



i=l 



flo + aiT, 



il~2^+Sa 



MlAp^^)- (22) 



The matrix elements of the operator ML^Li{pYi) = 
JL'{pri)[YLi{ri)ai]^ (with L' and cr coupled to L) are 
given in Appendix O 



B. Parity constraints 



The different multipoles in Eqs. (120])-(|22j) have well- 
defined parity 11, which can be deduced from the defi- 
nitions given in Appendix [Bl Eqs. (jB6p - (jB8[) . and the 
transformations under parity of 



n(v) 



-1, ii{Ylm) ^ {-i)\ n(Yiii)-(-i) 



and the parity of axial- vector one-body currents n( J"^) = 
-1-1. For elastic scattering, where the initial and final 
states of the nucleus are identical (J = Ji = Jf), only 
the multipoles with positive parity (H = +1) contribute 
to the structure factor, so that wc have 



n(/:i) = (-1 
U{Tt') = (-1 



,i+l 



,L + 1 



-1)^ 



L odd, 
L odd, 
L even . 



Hence, for elastic scattering only the odd-L multipoles 
of the longitudinal and transverse electric operators and 
only the even-L multipoles of the transverse magnetic 
operator contribute. This is also the case for inelastic 
scattering between initial and final states of the same 
parity. For inelastic scattering involving different parity 
states, the above constraints get reversed. 



C. Time-reversal constraints 

For elastic scattering, time-reversal invariance also 
constrains the multipoles that contribute to the struc- 
ture factor. We can write the reduced matrix elements 
of the sum over one-body operators OL{i) as [TjI 

A 

{J\\Y,OL{^)\\J) 



j,j' 



^jij,j'){{j\\OL\\j') + i-iy-'' (j'wolWj) 



(23) 



where ^j(j, j') denotes the one-body density matrix, and 
the sum is over single-particle total angular momenta 
J, j' (for simplicity, we have suppressed the sums over ra- 
dial quantum numbers n, n' and orbital angular momenta 
1,1'). Therefore, the symmetry properties of the matrix 
elements under exchange of initial and final states de- 
termine the allowed L contributions to elastic scattering. 
The relevant operator for SD WIMP-nucleus scattering is 
Ml,l'i whose matrix elements are given in Appendix [Cl 
They transform as 



(n'rij'||Mi,i(pr,)||n;ij) 



(_1)J+/ (nzij||Mi,i(pr, 



{n'l'\j'\ 



ML,L±l{pv^)\\nl-3) 



\n'l'\]') 



JvK 



(24) 



(25) 



Therefore, from Eq. ((23|) it follows that only the multi- 
poles with Ml^l±i contribute to elastic scattering. Con- 
sidering the different multipoles in Eqs. ()20|) -(|22 |) . we 
thus have 



i^-iy-i {nl-miL,L±i{pv^)\\n'l'-j') 



(Jlirr^^llj) 



0, 



(26) 
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FIG. 1. (Color online) Comparison of calculated spectra of 
^^®Xe and ^^^Xe with experiment. 



FIG. 2. (Color online) Comparison of the calculated ^^^I spec- 
trum with experiment. 



so that the transverse magnetic multipoles do not con- 
tribute to elastic scattering. 



D. Structure factor for elastic SD scattering 

As a resuh, the structure factor for elastic SD WIMP 
scattering off nuclei is given by [5| 

SAiP)= E (l('^ii^i(p)ii-^)l' + l('^iiTL"(p||'^)l'). 

L odd 

(27) 
and only odd-i longitudinal and electric transverse mul- 
tipoles contribute. 



IV. RESULTS 

A. Spectra 

The calculation of the structure factors requires a re- 
liable description of the nuclei involved in the scattering 
process. We perform state-of-the-art large-scale shell- 
model calculations of the nuclear states using the code 
ANTOINE [43 • For each nucleus, we solve the many- 
body problem in an appropriate valence space, which de- 
pends on the nuclear mass region. In all calculations, we 
use nuclear interactions that have been previously em- 
ployed in nuclear structure and decay studies. To test 
the quality of the structure calculations, we first com- 
pare the theoretical with the experimental spectra for all 
relevant isotopes. 



neutrons comprises the Ogj/2, 1^5/21 1'^3/2j 2si/2j and 
O/111/2 orbitals on top of a ^°°Sn core. For ^"^^Xe we per- 
form an exact diagonalization in this space. However, 
in order to make the calculations feasible for ^^^Xe, the 
number of particle excitations from the lower- lying 0177/2, 
1^5/2 orbitals into the 1^3/2, 2si/2, and Oft.11/2 orbitals 
was limited to three. With these restrictions the matrix 
dimension for this space is 3.5 x 10^. Similarly, for ^•^^I 
the number of excitations into the the 1^3/2, 2si/2, and 
Qhii/2 orbitals was limited to four, leading to a matrix 
dimension of 4.3 x 10^. For this valenc e sp ace we have 
used the so-called GCN5082 interaction [4^113, which is 
based on a G-matrix with empirical adjustments, mainly 
in the monopole part, to describe nuclei within this re- 
gion. The same interaction and valence space have been 
used to stu(iy nuclear structure and double-beta decays 
in Refs. ^^. 

Figure [1] shows the excitation energies of the lowest- 
lying states of ^^^Xe and -'^^-'^Xe in comparison with exper- 
iment (all energies are measured from the ground state). 
These spectra have been previously presented in Ref. [l6| . 
In Fig. [21 we show the spectrum of ^^^I. For all three 
cases, the experimental ground state and the overall or- 
dering of the excited states are very well described. This 
represents a clear improvement with respect to previ- 
ous work [2^, and validates the interaction and valence 
space used. Note that for ^^^I the spin and parity as- 
signment for some experimental states are not known. 
These states are absent in our calculated spectra, which 
suggests that they have significant contributions from or- 
bitals lying outside the valence space considered in the 
present calculations. 



'Xe, 



"■Xe, 



For the heaviest nuclei for SD WIMP scattering, ^^^Xe, 
-'^^^Xe and ^'^'^1, the valence space for both protons and 



'Ge 



For ''^Ge, the valence space for both protons and neu- 
trons comprises the 1^3/2, O/5/2, lpi/2, and 0^9/2 or- 
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FIG. 3. (Color online) Comparison of calculated '^^Ge spectra, FIG. 4. (Color online) Comparison of calculated spectra of 



using Int. 1 and Int. 2 interactions (for details see text), with 
experiment. 



bitals on top of a ^^Ni core. The calculations are per- 
formed in the complete space. We compare results for 
two different interactions, the so-called GCN2850 inter- 
action Wa|43| (Int. 1 in the following) and the RG inter- 
action [50] (Int. 2). Both are also based on a G-matrix, 
with mainly monopole empirical adjustments for this re- 
gion. They have been eniployed in beta and double-beta 
decay studies, Refs. [H, ll^ for Int. 1 and Refs. [13, [HI 
for Int. 2. The former was also used in a smaller valence 
space for the description of ''^Ge in Ref. [l3j . 

In Fig. |3] we compare the resulting spectra with ex- 
periment. We find that the ground state and the overall 
ordering of states is much better reproduced by the Int. 2 
interaction. In particular, the structure of three of the 
lowest-lying states and the gap between them and the 
higher-lying states are well described. In contrast, the 
Int. 1 interaction predicts a l/2~ ground state, in dis- 
agreement with experiment, and the general spacing of 
the spectrum is not well described. Consequently, the 
Int. 2 interaction will be the preferred one in this work. 
Nevertheless, we will also keep the Int. 1 case, in order 
to study the sensitivity of the structure factor to the dif- 
ferent nuclear interactions. It is important to note that 
the first exited state, which is a 5/2+ state, is at too 
high excitation energy in both calculations. This sug- 
gests that an extended valence space, probably including 
the higher-lying 1^5/2 orbital, is needed to account for 
this state. This was also observed in Ref. [13]. A reliable 
description of the 5/2+ state will be crucial for the study 
of inelastic scattering off ^'^Ge. 
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FIG. 5. (Color online) Comparison of calculated spectra of 
^^Al and ^^Si with experiment. 



tions in this valence space are easily performed. In pre- 
vious works [ll[I3,[ll,[2S[ill, the USD interaction [52 1 
was employed. This interaction consists of a best fit to 
selected nuclei in this mass region. Here, we use the more 
recent USDB interaction [53|, which is an improved ver- 
sion of USD. The difference between the two interactions 
is small, see Sect. IIVBI In Figs. U] and [5] the positive- 
parity excited states of all four nuclei are shown com- 
pared to experiment (in the sd shell only positive parity 
states can be obtained). The agreement with experiment 
is very good in all cases, both for the ordering and the 
quantitative reproduction of the excitation energies. 



'F, ^'-"Na ,^^Al, ^""Si 



B. Spin expectation values 



The valence space of the four lighter nuclei ^^F, ^'^Na, 
^^Al, and ^^Si is the sd shell, which comprises the 0^5/2, 
lsi/9, and Od-i/2 orbitals, with a ^^O core. Full calcula- 
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In the limit of low momentum transfer, p = 0, the 
structure factor for elastic SD WIMP scattering is given 
by the proton and neutron spins Sp = X]i=i ""1/2 and 
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TABLE III. Calculated spin expectation values for protons (Spjand neutrons (S„> of ^^^'"^Xe, ^^^I, ''^Ge, ^^Si, ^''Al, ^^Na, 
and ^^F, compared to the previous calculations of Refs. [la. ll7H23l | . 



S„ = J2i=i ""1/2 in the nucleus Q: 



(2J+1)(J+1) |^ ^'V«\^^ 'V« \|2 

— (ao-hai)(Sp) 4- (ao-aJ(S„) , 

4:TrJ 



Sa{0) = :i^|(ao + a[){J\\Sp\\J) + {ao - a[){J\\S,,\\J)\^ 
471- 

(28) 

(29) 

where a'l = ai(l -|- Sai) includes the effects from chiral 
2b currents. The spin expectation values are defined as 
(S„,p) = (JM = J|S3_p|JM = J). 

We list our calculated spin expectation values (S„^p) 
in Table IIIII in comparison to previous calculations. As 
expected for odd-mass nuclei with even number of pro- 
tons (i2943ixe, 73Ge, and ^ssi) |(S„)| > |(Sp)|, while for 
odd-mass nuclei with an even number of neutrons (^^F, 
23Na ,27A1, and ^^n) |(S„)| < |(Sp)|. As a result, the 
WIMP coupling to the even species will be suppressed. 
Moreover, the sensitivity to the precise value of the even 
species spin is very weak when chiral 2b currents are in- 
cluded. This is shown in Sect. IIVCI Chiral 2b currents 
lead to an interaction of neutrons and protons that over- 
whelms the direct WIMP coupling to the suppressed spin 
expectation value, so that the structure factors are al- 
most entirely determined by the dominant (S„/p) (for 
odd neutron/proton isotopes). 

The spin expectation values of the lighter nuclei, ^^F, 
^^Na, ^^Al, and ^^Si in Table Hill are very close to those of 
Refs. [ll,[i3,[ll,[2l[2l| due to the similarity of the USD 
and USDB interactions. This indicates that the struc- 
ture for these nuclei is under good control. For ^■^Ge we 
find a weak sensitivity of the dominant (S„) value com- 
paring the preferred Int. 2 interaction ("This work") to 
the Int. 1 interaction. This range is smaller than the one 
in previous calculations of Refs. [13, [l3, [l^ [221 ; suggest- 
ing that the latter may have an even larger variation in 
the spectra due to truncations or deficiencies in the inter- 
actions used. Also for the heavier nuclei, ^^^'^'^^Xe, and 



^'^'^I, we have performed calculations in the largest spaces 
to date and with tested interactions. For ^^^'^"^^Xe, the 
comparison to previous results is discussed in detail in 
Ref. [3. For the dominant (S„) values for I29,i3ixe, 
and the dominant (Sp) value for ^^^I, the difference to 
previous calculations of Refs. [ij, [20, [22, [23 is about 
25% (and 55% for ^^^Xe). We attribute these differences 
to the sizable truncations of the valence spaces in those 
calculations and because the interactions used have not 
been as well tested. 



C. Structure factors 

1. Isoscalar/isovector versus proton/neutron 

The structure factor Sa(j>) can be decomposed in 
terms of its isoscalar and isovector parts Sij{p), charac- 
terized by the isoscalar and isovector couplings aq and ai : 

Sa{p) = al Sooip) + aoaiSoiip) + a\ S'ii(p) . (30) 

However, it is common in the literature to use the struc- 
ture factors Spip) and Snip), which are referred to as 
"proton-only" and "neutron-only" , respectively. They 
are defined by the couplings ao — ai — 1 ("proton-only") 
and ao = — ai = 1 ("neutron-only") and are thus related 
to the isoscalar and isovector structure factors by 

Sp{p) = Sooip) + Soiip) + Snip) , (31) 

Snip) = Sooip) - Soiip) + Snip) . (32) 

The origin of the "proton/ncutron-only" structure fac- 
tors can be understood from Eq. ([^5]) . When 2b cur- 
rents are neglected, at p = the "proton/neutron-only" 
structure factors are determined entirely by the pro- 
ton/neutron spin expectation values. Moreover, when 
the higher-order isovector parts in lb currents are ne- 
glected, this separation also holds for p > 0. Because 
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FIG. 6. (Color online) Structure factors Sp{u) (solid lines) 
and Sn{u) (dashed) for ^'^^Xe (top panel) and ^^'^Xe (bottom 
panel) as a function of m = p^b^/2. The harmonic-oscillator 
lengths are b = 2.2853 fm and b = 2.2905 fm for ^^^Xe and 
^ ^Xe, respectively. Results are shown at the lb current level, 
and also including 2b currents. The estimated theoretical un- 
certainty is given by the red {Sp{u)) and blue {Sn{u)) bands. 



for odd-mass nuclei there is a clear hierarchy of the 
spin expectation values (with either |(S„)| ^ |(Sp)| or 
|(Sp)| >• |(S„)|), the proton/neutron decomposition is 
useful to capture the dominant parts of Sa^p)- For this 
reason, and because it is common experimentally, we will 
also largely consider the proton/neutron decomposition 
here. This is merely a convenient choice of ao,ai cou- 
plings, but the notation "proton/neutron-only" is mis- 
leading, because it does not imply that the coupling is 
to protons/neutrons only. Strong interactions between 
nucleons in 2b currents, as well as the isovector nature 
of pseudo-scalar and other Q^ lb currents, mean that 
WIMPs effectively couple to protons and neutrons in nu- 
clei. In fact, with 2b currents, both Sp{p) and Sn{p) are 
determined by the spin distribution of the odd species. 

In the following, we present structure factors as a 
function of m = p^b'^/2 with harmonic-oscillator length 
b = {h/mujY^^ and huj = (45^-1/3 - 25^-2/3) MeV. 
When 2b currents are included, we provide theoretical 
error bands due to the uncertainties in WIMP currents 
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FIG. 7. (Color online) Decomposition of the isovector struc- 
ture factor Sii(u) for ^^^Xe. At the lb current level, the 
full result (solid, black lines) and the contributions from 
transverse electric (dashed, blue) and from longitudinal (dot- 
dashed, green) multipoles are shown. The top panel gives 
also the full lb plus 2b current result (red band), while the 
middle/bottom panels show the lb plus 2b results when only 
transverse/longitudinal multipoles are included (blue/green 
band). The bands give the estimated 2b-current uncertainty. 



in nuclei, see Table |TT1 This takes into account the un- 
certainties in the low-energy couplings C3 , C4 and in the 
density range p = 0.10...0.12fm~ . 

For ^^^Xe and ^^iXe the predicted isoscalar/isovector 
structure factors S'oo(u), 5'qi(u), and S'ii(u) were dis- 
cussed in detail in Ref. |16|, and they were compared 
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FIG. 8. (Color online) Structure factors Sn{u) (top panel) 
and Sp{u) (bottom panel) for ^'^^Xe. At the lb current level, 
the full result (solid, black lines) are compared with the con- 
tributions from transverse electric (dashed, blue) and from 
longitudinal (dot-dashed, green) multipoles. 



to the previous calculations of Refs. [20, l2j{ (see also 
Sect. lIVB]) . Here, we present in Fig.[S]tlie proton/neutron 
structure factors Sp{u). At the lb current level, the re- 
sults at p = are determined by the spin expectation 
values. Chiral 2b currents provide important contribu- 
tions to the structure factors, especially for p < 100 MeV, 
where we find in Fig. [5] a significant increase of Sp{u). 
This is because with 2b currents, neutrons can contribute 
to the "proton-only" (oq = ai = 1) coupling due to the 
axial 6ai contribution in Eq. (|29| . For Sn{u), 2b cur- 
rents lead to a small reduction in the structure factor, 
depending on the momentum transfer. This is caused 
by the combined effect of the axial 6ai and the pseudo- 
scalar daf[p) contributions. To better understand how 
these different contributions enter, we study a multipole 
decomposition of the structure factors. 



2. Multipole decomposition 

In Fig. [7] we show the transverse/longitudinal decom- 
position of the results with lb as well as lb plus 2b cur- 
rents for the isovector structure factor 5*11 (m) of ^^^Xe 
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FIG. 9. (Color online) Decomposition at the lb current level 
of the ^^^Xe structure factor Sn{'u) (solid, black line) in L = 1 
(dashed, violet) and L — 3 (dot-dashed, orange) multipoles. 



(the long-range 2b currents are isovector). The different 
2b current contributions can be clearly seen in Fig. [T] In 
the middle panel, where only the transverse electric mul- 
tipoles are taken into account, 2b currents reduce the 
lb result due to the negative axial Sai values. Because 
(5a 1 is momentum independent, the relative reduction is 
independent of u. The bottom panel shows the longitu- 
dinal multipoles, where both axial Sai and pseudo-scalar 
(5af (p) 2b current contributions enter. At low momen- 
tum transfers we find a reduction of the structure factor, 
driven by 6ai , which turns into an enhancement at higher 
momentum transfer due to (5af (p). Both contributions 
are comparable at u ^ 3 corresponding to p ~ 200 MeV. 
In the upper panel, the full lb plus 2b band is given, 
where the final reduction or enhancement over the lb re- 
sult, for a given u value, depends on the relative impact 
of the transverse electric and longitudinal multipoles. 

It is interesting to study the transverse/longitudinal 
decomposition at the lb level, as shown in Fig. |S] for 
Sn{u) and Sp{u) of ^■^^Xe. While both multipoles con- 
tribute to Sn(u) (their relative importance depends on u), 
Sp{u) is completely dominated by the longitudinal mul- 
tipoles except at p = 0. In ^"^^Xe almost all of the spin 
is carried by neutroiEis, so >S'p(0) is very small at the lb 
level. However, for p > the (isovector) pseudo-scalar 
currents allow neutrons to contribute to Sp{u), leading 
to a steep increase in the longitudinal contribution to 
Sp{u). Because pseudo-scalar currents only contribute 
to the longitudinal multipoles, the transverse part from 
the protons also remains very small for p > 0. 

Another way to decompose the structure factors is in 
terms of the different L values of the multipoles. Be- 
cause the ground state of ^^^Xe is 1/2+ , only L — 1 con- 
tributes. For ^^^Xe, with a 3/2"*" ground state, L = 1 and 
L — 3 multipoles enter (even-L lEnultipoles are forbidden 
due to parity, see Sect. IIIIBp . The L decomposition of 
the ^■^^Xe structure factor S'„(u) is shown in Fig. [51 for 
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FIG. 10. (Color online) Structure factors Sp{u) (solid lines) 
and S„{u) (dashed) for '^^Ge calculated using the Int. 1 
(GCN5028, top panel) and the Int. 2 interaction (RG, bot- 
tom panel) as a function of u = p^b^/2 with b = 2.1058 fm. 
Results are shown at the lb current level, and also including 
2b currents. The estimated theoretical uncertainty is given 
by the red {Sp{u)) and blue {Sn{u)) bands. 



simplicity at the lb current level. We observe that the 
L = 3 multipoles dominate for 1.5 < m < 5. As a result, 
the structure factors fall off considerably more slowly for 
^^^Xe compared to ^^^Xe, where only L = 1 contributes. 
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FIG. 11. (Color online) Structure factors Sp{u) (solid lines) 
and Sn{u) (dashed) for ^ I as a function oi u — p b /2 with 
b = 2.2801 fm. Results are shown at the lb current level, and 
also including 2b currents. The estimated theoretical uncer- 
tainty is given by the red {Sp{u)) and blue (5,1(11)) bands. 



I 



127 j^ 19p^ 23^^^ 27^;_ 29^^ 



In Figs, [m [121 and [131 we show the structure fac- 
tors Sn{u) and Sp{u) for ^^"^l, ^^F, ^SNa, ^^Al, and ^SSi 
at the lb current level and including 2b currents. The 
dominant structure factor is the one for the odd species. 
Therefore, for ^^Si Sn{u) dominates, while for the other 
isotopes Sp{u) is the main component. All the features 
discussed for ^"^^Xe in Sect. II V C2| translate to these iso- 
topes as well: The structure factors for the non-dominant 
"proton/neutron-only" couplings is strongly increased 
when 2b currents are included. For the dominant struc- 
ture factor, 2b currents produce a reduction, by about 
10 — 30% at low momentum transfers, which at large u 
can turn into a weak enhancement due to the 2b current 
contribution to the pseudo-scalar currents. This is most 
clearly seen for ^^F in the top panel of Fig. [1^ where we 
also show the isoscalar/isovector structure factors Sqo (u) , 
S'oi(u), and Sii{u). Note that the structure factor S'oi(m) 
vanishes at the point where Sp{u) and S'n(w) cross. 



3. 



'Ge 



Figure [TOl shows the structure factors for ^^Ge for the 
different Int. 1 and Int. 2 interactions (the latter is pre- 
ferred based on the spectra, see Fig. ^ . The structure 
factor Sn{u) differs by less than 10% between the two 
interactions. At the lb current level, Sp{u) for low mo- 
mentum transfers is substantially smaller for Int. 1, due 
to the very small (Sp) value. However, when 2b cur- 
rents are included, also for Sp{u) the contributions from 
neutrons are dominant, which translates to similar struc- 
ture factors for the two interactions. This is because of 
the similar (S„) values (see Table HID) combined with the 
neutron-proton coupling through 2b currents. 



V. CONCLUSIONS AND OUTLOOK 

This work presents a comprehensive derivation of SD 
WIMP scattering off nuclei based on chiral EFT, includ- 
ing one-body currents to order Q^ and the long-range 
Q'^ two-body currents due to pion exchange, which are 
predicted in chiral EFT. Two-body currents are the lead- 
ing corrections to the couplings of WIMPs to single nu- 
cleons, assumed in all previous studies. Combined with 
detailed Appendices, we have presented the general for- 
malism necessary to describe both elastic and inelastic 
WIMP-nucleus scattering. 
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FIG. 12. (Color online) Structure factors for ^^F as a 
function oi u = p^h^ /2 with b = 1.7608 fm. Top panel: 
isoscalar/isovector Soo(u) (solid line), S'oi(u) (dashed), and 
S'ii(it) (dot-dashed) decomposition. Bottom panel: pro- 
ton/neutron Sp{u) (solid line) and Sn{u) (dashed) decom- 
position. In both panels results are shown at the lb current 
level, and also including 2b currents. The estimated theoret- 
ical uncertainty is given by the red (Sii(u), Sp{u)) and blue 
(S'oi(ii), S„{u)) bands. 




0.001 b- 
0.1 



3 



0.001 r 



0.0001, 



S (u) lb currents 



S (u) lb currents 

S (u) lb + 2b currents 

p 

S (u) lb + 2b currents 








We have performed state-of-the-art large-scale shell- 
model calculations for all non-zero-spin nuclei relevant 
to direct dark matter detection, using the largest valence 
spaces accessible with nuclear interactions that have been 
tested in nuclear structure and decay studies. The com- 
parison of theoretical and experimental spectra demon- 
strate a good description of these isotopes. We have cal- 
culated the structure factors for elastic SD WIMP scat- 
tering for all cases using chiral EFT currents, including 
theoretical error bands due to the nuclear uncertainties 
of WIMP currents in nuclei. Fits for the structure factors 
are available upon request. 

We have studied in detail the role of two-body currents, 
the contributions of different multipole operators, and 
the issue of proton/neutron versus isoscalar/isovector de- 
compositions of the structure factors. The long-range 
two-body currents reduce the isovector parts of the struc- 
ture factor at low momentum transfer, while they can 
lead to a weak enhancement at higher momentum trans- 



FIG. 13. (Color online) Structure factors Sp{u) (solid lines) 
and Sn{u) (dashed) for ^'^Na (top panel), ^^Al (middle panel), 
and Si (bottom panel) as a function of u = p b /2, 
with harmonic-oscillator lengths b — 1.8032 fm (^'^Na), b = 
1.8405 fm (^'^Al), and b = 1.8575 fm (^^Si). Results are shown 
at the lb current level, and also including 2b currents. The 
estimated theoretical uncertainty is given by the red (Sp{u)) 
and blue {Sn{u)) bands. 



fers. Moreover, we have shown that for odd-neutron 
(odd-proton) nuclei, two-body currents lead to a sig- 
nificant increase of the "proton-only" ("neutron-only") 
structure factors, because of strong interactions between 
nucleons through two-body currents that allow the odd 
species carrying most of the spin to contribute. This im- 
plies that WIMPs effectively couple to protons and neu- 
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trons in nuclei, so that the notation "proton/neutron- 
only" is misleading. In fact, with 2b currents, both 
"proton/neutron-only" structure factors are determined 
by the spin distribution of the odd species. 

Future improvements of the nuclear physics of dark 
matter detection includes developing shell-model inter- 
actions based on chiral EFT, where the present frontier 
are semi- magic nuclei up to the calcium region J54-58J, 
ab-initio benchmarks for the lightest isotope ^^F, and ex- 
panding the valences spaces (especially for germanium). 
In addition, a full treatment of the one- and two-body 
currents would require to renormalize them to the va- 
lence space of the many-body calculation, which can lead 
to additional contributions to the currents. This and go- 
ing beyond the normal-ordering approximation will be 
pursued in future work. Moreover, we plan to investi- 
gate other responses [13] based on the same large-scale 
nuclear-structure calculations presented here. 
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Appendix A: 
Calculation of the effective one-body current Jf^b 



We calculate the normal-ordered lb part of 2b currents 
by summing the second nucleon over occupied states of a 
spin and isospin symmetric reference state or core, which 
we take as a Fermi gas: 



terms of the C3 term in Eq. 



Ji,^r(c4term) 



F2 2 
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C4 



-V 



Am J (27r)3 
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where /f (p, P) , /f (p, P) , and I^ (p, P) are given by the 
following integrals 



cff 
i,2b 



J 



(1-P,,)J 



(Al) 



where the sum is over occupied states, J?^ is the 2b cur- 
rent defined in Eq. ([5]), and P^ is the exchange operator. 
In this approximation, the momentum k in the direct (d) 
and exchange (ex) contributions is given by 



^(P- 



k- = /^';.kd = p,-9^. 



Pj + PjJ 
p 
2 



p 

2 ' 



(A2) 
(A3) 



where we have used that the initial and final momenta of 
the nucleon in the occupied state are identical, p^ — p'-. 

The non-vanishing contributions to J^^^ can be 
grouped into three terms, arising from the first line of 
Eq. ([5]) (c4 term) as well as from the direct and exchange 
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Appendix B: 
Derivation of the structure factor Sa{p) 

We start from the Lagrangian density for spin- 
dependent WIMP-nucleus scattering Eq. ([!]). WIMPs 
are expected to be nonrelativistic with velocities of the 
order v/c ^ 10^'^, so the time components of the currents 
can be neglected. Evaluating the Lagrangian density be- 
tween initial and final states leads to 

(/|£f 10 ^-^Id^re-P-X/TTSX.J^.W, (Bl) 

where e~*P'""x/77^Xi — (x/l j(r) \xi} represents the ma- 
trix element of the leptonic current of the WIMP and 
j4j(r) that of the hadronic current. 

We can expand the leptonic current in terms of spher- 
ical unit vectors 15911: 



X/77 Xi 



=>i/- e-'P"" — 1 o-*P'' 



Ig-^pr^ Y^ l^ele-'P", (B2) 

A=0,±1 

with spherical unit vectors with z-axis in direction of p 

e±i = T^=(epi ±iep2) Gq = -p-, , (B3) 

v2 IPI 
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'±1 = T^('i ± «^2) lx=n = k' 



(B4) 



We can also expand the product e^ e *p'' in Eq. (|B2[) in 
a multipole expansion |59| . This leads to 

G 



-^ (J/Af/I f ^ VM2L+l){-^)%Clo{p) 

J2 v/M2rTT)(-z)^ 



L>1 



E l>^[rex{p) + ^Trfip)]]\J:M.) , (B5) 



A=±l 



where \JiMi), \JfMf) denote the initial and final states of 
the nucleus, p = |p|, and the electric longitudinal, electric 
transverse, and magnetic transverse multipole operators 
are defined by [SQlj 



P 



V[jLipr)YLM{nr)] -J^W, 
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71T/'(P) ^Jd'r [jUpr)YtUnr)] ■ J^(r) , (B8) 

with spherical Bessel function jl (pr) . The vector spher- 
ical harmonics are given by 

YfLa(^r) = E (i'™lA|i'liM) YL',n{nr)ex . (B9) 



TnA 



Since J (r) = J2i=i Jf (i")'^(i" ~ i"j)j the multipole opera- 
tors can be written as a sum of one-body operators: 
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The structure factor Sa{p) is obtained from 



(/l>C 



SD 



by summing over the final neutrahno 



spin and over the nucleus final-state angular momentum 
projections, and by averaging over the initial configu- 
rations. It is thus useful to work with reduced matrix 
elements that do not depend on projection numbers: 



{JfMf\OLM\J^M,) 
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with 3j coefficients and where O is a tensor operator of 
rank L. This gives for the sum and average [59j 
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where we have assumed that the neutrahno spin is 1/2, 
and the cross terms vanish due to the orthogonalization 
properties of the 3j coefficients. For the sum over neu- 
trahno spin projections one has for ji^v = 1,2,3 

Si,Sf 

= E {x?{pf)xJ{Pf){r7'u4Hp^m{p^)h'ln,s), 

= Tri-ff'j^j'^j'') = -'iS^"' , (B15) 

which follows from the completeness relation 
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valid for nonrelativistic WIMPs. Combined, this gives 



the final result: 

^E i:i(/i4°i.>i 
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The specific form of the multipoles depends on the 
form of the WIMP currents Jf (r) . They contain either 
axial- vector terms [(Ti] or pseudo-scalar ones [(p • (Ti)p]. 
For axial- vector currents, the response will be propor- 
tional to the following operator 

Mt'Mpr^) = jL'(P^.)Yili(rO • a, , 

= JL'ipn) Y, {L'mlX\L'lLM) Yl'^{y,) a}^ , 



iL'{pri)\YL'{x 



(B18) 



Pseudo-Scalar currents, which are proportional to the 
momentum transfer p, only contribute to the longitu- 
dinal multipoles (see Eq. ([B5|). Moreover, in these we 
can replace (p • cri)p by p~<Ti, because of 



(p • (Ti)l> = p^CTi + p X (p X (Ti) , 



(B19) 



and the second term is perpendicular to p, so it vanishes 
for the longitudinal multipoles. As a result, pseudo-scalar 
currents can also be expressed in terms of Mffj^,{pVi). 

In summary, including chiral 2b currents at the 
normal-ordered one-body level in Eq. (J17p . we have for 
the multipoles 
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Note that the p'^/A\ terms cancel in the longitudinal re- 
sponse and only contribute to the transverse multipoles. 
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Appendix C: 
Reduced matrix elements of ML,L'{pvi) 

To calculate the structure factor, we need the ma- 
trix elements of the one-body operator M^^L'ip^i) = 
JL'{p'r'i)[YL'{ri) <Ti\^ between the single-particle states of 
the many-body basis used for the description of the nu- 
clear states. The reduced matrix elements can be ob- 
tained as a function of 3j and 9j symbols and matrix 
elements of the spherical Bessel functions j^' , 
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